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ABSTRACT 

We present the results of a photometric and spectroscopic study of the white dwarf can- 
didate members of the intermediate age open clusters NGC3532 and NGC2287. Of the nine 
objects investigated, it is determined that six are probable members of the clusters, four in 
NGC3532 and two in NGC2287. For these six white dwarfs we use our estimates of their 
cooling times together with the cluster ages to constrain the lifetimes and masses of their pro- 
genitor stars. We examine the location of these objects in initial mass-final mass space and 
find that they now provide no evidence for substantial scatter in initial mass-final mass relation 
as suggested by previous investigations. Instead, we demonstrate that, when combined with 
current data from other solar metalicity open clusters and the Sirius binary system, they hint at 
an IFMR that is steeper in the initial mass range 3MQ^Minit^4M0 than at progenitor masses 
immediately lower and higher than this. This form is generally consistent with the predictions 
of stellar evolutionary models and can aid population synthesis models in reproducing the 
relatively sharp drop observed at the high mass end of the main peak in the mass distribution 
of white dwarfs. 

Key words: stars: white dwarfs; galaxy: open clusters and associations: NGC2287; 
NGC3532 



1 INTRODUCTION 

Standard stellar evolutionary theory predicts the existence of a 
positive correlation between the main sequence mass of a star 
(M;$1OM0) and the mass of the white dwarf remnant left behind 
after it has expired. This is frequently refered to as the initial mass- 
final mass relation (IFMR). A secure and detailed knowledge of the 
form of the IFMR is important to a number of very active areas of 
astrophysical research. For example, the relation is a key ingredient 
of models of the chemical evolution of the Galaxy as it provides an 
estimate of the amount of gas, enriched with C, N and other met- 
als, a low or intermediate mass star returns to the ISM (e.g. Carigi, 
Colin & Peimbert 1999). Moreover, understanding the form of the 
IFMR is crucial to deciphering information locked up in the white 
dwarf luminosity functions of stellar populations (e.g. leffery et 
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al. 2007, Oswalt et al. 1996). The shape of the upper end of the 
IFMR is of special interest as it places limits on the minimum mass 
of a star that will experience a Type II supernova explosion. With 
robust constraints on this mass, for example, the observed diffuse 
SNe neutrino background can serve as an empirical normalisation 
check on estimates of the star formation history of the Universe 
(e.g. Hopkins et al. 2006). 

Unfortunately, due to the complexity of the physical processes 
occuring within stars, especially during the final stages of the stellar 
lifecycle, the form of the IFMR is rather difficult to ascertain from 
first principles. For example, the mass of the stellar core which ul- 
timately becomes the white dwarf is probably modified with each 
thermal pulse cycle on the asymptotic giant branch (AGB). The 
final remnant mass predicted by evolutionary models is therefore 
dependent on the length of time a star is assumed to spend on the 
AGB. However, evolution on the AGB is terminated by the removal 
of the stellar envelope so the predicted duration of this phase is sus- 
ceptible to the assumptions made about the rate at which envelope 
mass is lost (e.g. Iben & Renzini 1983). While significant inroads 
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are being made in tlie tlieoretical understanding of mass loss on tiie 
AGB (e.g. Bowen 1988, Wachter et al. 2002), a comprehensive and 
robust physical treatment remains elusive. 

Therefore, empirical data currently play a crucial role in im- 
proving our understanding of the form of the IFMR. Arguably the 
best way at present to obtain observational based constraints on the 
IFMR is through the study of the white dwarf members of open star 
clusters. The ages of these coeval populations can be determined 
from the location of the main-sequence turn-off or the lithium de- 
pletion boundary in the sequence of the lowest mass members (e.g. 
Pleiades; Stauffer et al. 1998). The progenitor star Ufetimes and ul- 
timately their masses can then be estimated by calculating the dif- 
ference between the age of the parent population and the cooling 
times of white dwarf members. It is worth pointing out, however, 
that as theoretical input is required at this latter stage, in the form 
of stellar evolutionary models, constraints acquired in this way 
are more appropriately termed semi-empirical (e.g. Weidemann & 
Koester 1983). 

With greater access to 8m class telescopes and the availabil- 
ity of improved instrumentation, the last few years has seen a 
flurry of new studies of the white dwarf candidate members of 
open star clusters aimed at constraining the form of the IFMR. 
For example, Claver et al. (2001) have presented the results of a 
detailed pholoniclric and spectroscopic study of five white dwarf 
members of Praesepe. Williams, Bolte & Koester (2004) have un- 
dertaken a spectroscopic study of the degenerate candidate mem- 
bers of NGC2168. Kalirai et al. (2005) have obtained multi-object 
spectroscopy of well over a dozen white dwarf candidate mem- 
bers of NGC2099. Dobbie et al. (2004,2006a) and Casewell et al. 
(2009) have identified and spectroscopically analysed a number of 
new white dwarfs in Praesepe. Williams & Bolte (2007) have in- 
vestigated degenerate candidate members of both NGC6633 and 
NGC7063, while Kalirai et al. (2008) present spectroscopy of white 
dwarf candidate members of two older open clusters, NGC7789 
and NGC6819. Most recently, Rubin et al. (2008) have unearthed 
six white dwarf candidate members of the intermediate aged clus- 
ter NGC1039. Combined, these studies have increased the number 
of data points on the semi-empirical IFMR by over a factor of two. 
However, only the work of Williams, Bolte & Koester (2004, 2008) 
and Rubin et al. (2008) has provided a substantial number of new 
points in the Minit^4M0 regime. The upper portion of the IFMR 
remains rather poorly constrained, with large scatter amongst the 
data, the bulk of which currently comes from only three clusters, 
NGC1039, NGC2516 and NGC2168. Indeed, to further compli- 
cate the situation, the metalicity of NGC2168 is subsolar by ap- 
proximately a factor two, and thus the data points from this cluster 
might not be expected to occupy the same region of initial mass- 
final mass space as those from solar metalicity clusters such as the 
Pleiades, NGC2516 and NGC1039 (e.g. Marigo 2001). 

In a bid to improve the current state-of-affairs, we have re- 
cently focused our efforts on open clusters with r;$300Myrs, cor- 
responding to the lifetime of a M«4Mo star, and with metalici- 
ties close to solar so that the influence of this parameter on the 
data from cluster to cluster is minimal or even negligible. Here we 
present the results of a new investigation of the nine white dwarf 
candidate members of NGC2287 (M41) and NGC3532. In the next 
section we review the parameters of these two clusters and the work 
previously undertaken on their white dwarf populations. Next, we 
describe the acquisition and analysis of new photometry and spec- 
troscopy for these objects and then use these to re-assess cluster 
membership status. FinaUy, we examine those white dwarfs with 



parameters consistent with being members of the two clusters in 
tiie context of tiie IFMR. 

2 TWO OPEN CLUSTERS RIPE FOR A NEW 8M STUDY 

2.1 White dwarf candidate members of NGC3532 

NGC3532 (ll''05'" -58°45', J2000.0) is a comparatively rich and 
nearby open cluster yet has been relatively little studied, partic- 
ularly in the last two decades. This is probably because it ap- 
pears in projection against the Galactic plane (b~1.3°). Despite 
this low Galactic latitude, extinction along the line of sight to 
the cluster is low. For example, Fernandez & Salgado (1980) and 
Meynet, Mermilliod & Maeder (1993) use UBV photometry to 
determine E(_B — V^)~0.04, from Stromgren data Eggen (1981) 
estimate E(b - ;y)=0.023 (which equates to E(B - F)«0.03) and 
Nicolet ( 1 98 1 ) find E(B - V)=0.052±0.0 1 using Geneva photom- 
etry. Moreover, NGC3532 appears to have a metalicity close to the 
solar value. From a photometric investigation of the giant members, 
Claria & Lapasset (1988) found no evidence for a significant UV 
excess or cyanogen anomaly while Twarog et al. (1997) concluded 
from an independent DDO based study that [Fe/H]=-0.02. 

While it is unfortunate that there are less than a handful of 
age estimates for NGC3532 in the recent literature, at least the few 
that are available appear to be consistent. For example, Meynet, 
Mermilliod &. Maeder (1993) compared isochrones generated from 
the stellar models of SchaUer et al. (1992), which include a mod- 
erate level of convective core overshooting, and UBV photometry 
to obtain t=316M5TS. Koester & Reimers (1993), using moderate 
core overshoot models, determined an age of r=302±154Myrs for 
a turn-off absolute magnitude of Mv=-0.75±0.25. Kharchenko et 
al. (2005) estimate an age of r=282Myrs based on photometry from 
the 2MASS Point Source Catalogue (PSC; Skrutskie et al. 2006) 
and the stellar modelling of the Padova group (e.g. Girardi et al. 
2000). Therefore, it appears reasonable to conclude that the age of 
NGC3532 lies within the range T=300±25Myrs. 

Robichon et al. (1999) used Hipparcos data to determine the 
cluster distance modulus to be m-M=8.04lQ 32 . This is notably 
less than the recent estimate of m-M=8.61 obtained from 2MASS 
PSC data (Kharchenko et al. 2005). Nevertheless, these two de- 
terminations bracket the results from most other studies. Fernan- 
dez & Salgado (1980) estimate m-M=8.45±0.27, Eggen (1981) 
find m-M=8.5±0.25 while Meynet, Mermilliod & Maeder (1993) 
conclude that m-M=8.35. Reimers & Koester (1989) exploited this 
relative proximity and the low levels of foreground reddening, us- 
ing deep UV and red Schmidt plates to identify seven blue can- 
didate white dwarf cluster members. Low resolution spectroscopy 
obtained with the 3.6m ESO telescope and the Faint Object Spec- 
trograph and Camera (EFOSC) confirmed the degenerate nature of 
three of these stars. A subsequent extended photometric and spec- 
troscopic survey of NGC3532 unearthed a further three white dwarf 
candidate members (Koester & Reimers 1993). 

While analysis of the EFOSC spectroscopy confirmed that 
these white dwarfs are young (Tcfj ~20000-30000K), as would be 
expected if they are members of NGC3532, due to the low signal- 
to-noise of this data blueward of 4000A, the uncertainties in deter- 
minations of their surface gravities and hence in estimates of their 
masses and cooling times, remain large (e.g. see Koester & Reimers 
1993). The location of these white dwarfs in initial mass-final mass 
space as revealed by a number of recent studies e.g. Ferrario et al. 
(2005), Catalan et al. (2008) is suggestive of significant dispersion 
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in the IFMR. A large intrinsic scatter in tlie relation would have 
important implications for our understanding of stellar evolution 
and pose a huge complication for investigations which rely on the 
IFMR being close to monotonic (e.g. Jeffery et al. 2007). However, 
before firm conclusions are drawn on this matter, better photomet- 
ric and spectroscopic data should be acquired for these stars so that 
the cluster membership status of each can be soundly examined and 
more stringent limits placed on their masses and cooUng times. 



2.2 White dwarfs candidate members of NGC2287 (M41) 

An early detailed investigation of the moderately populated 
NGC2287 (06''46™ -20°45', J2000.0) determined the distance to 
the cluster to be D=725pc and levels of extinction along this line of 
sight to be very low, E{B - V')=0.01 (Cox 1954). These estimates 
are corroborated by the results of more recent studies where state- 
of-the-art theoretical isochrones have been fit to the observed clus- 
ter sequence. For example, based on an analysis of UBV photom- 
etry of cluster members Meynet, Mermilliod & Maeder (1993) de- 
rive a distance modulus of m-M=9. 15 and estimate E(S — V)=OM. 
Kharchenko et al. (2005) determine m-M=9.30 and E(B-V)=0.03 



based on a study of near-lR data from the 2MASS PSC (Skrutskie 
et al. 2006). Sharma et al. (2006) estimate m-M=9.26 and E(B- 
V)=0.01 from BVI imaging obtained with the Kiso Schmidt tele- 
scope. 

Recent age determinations for the cluster, which employ stel- 
lar models including a moderate level of convective core over- 
shoot, are at a reasonable level of agreement with one and other. 
For example, Meynet, Mermilliod & Maeder (1993) conclude that 
the cluster is T=240Myrs, while Harris et al. (1993) determine 
an age of T=200Myrs. Kharchenko et al. (2005) estimate the 
age of NGC2287 to be r=280Myrs, while Sharma et al. (2006) 
favour r=250Myrs. Thus the weight of evidence favours the age of 
NGC2287 to he within the range r= 243±40Myrs. However, met- 
alicity estimates for NGC2287 are somewhat scattered. In their cat- 
alogue of spectroscopic stellar abundance determinations, Cayrel 
de Strobel et al. (1992) quote values of [Fe/H]=0.00 and [Fe/H]=- 
0.25 for the cluster members HD49091 and HD49068 respectively. 
Cameron (1985) used UBV photometry of 14 members to derive 
[Fe/H]=-i-0.065 while from Stromgren photometry of F star mem- 
bers, Nissen (1988) determined [Fe/H]=-0.10. Collectively these 
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Table 1. Details of the nine white dwarf candidate members of NGC3532 and NGC2287. Masses and cooling times for each 
star have been estimated using the mixed CO core composition "thick H-layer" evolutionary calculations of the Montreal Group 
(e.g. Fontaine, Brassard & Bergeron 2001). The eiTors in absolute magnitudes, masses and cooling times shown here are derived 
by propagating more realistic uncertainties in effective temperature and surface gravity of 2.3% and O.OVdex respectively. 
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results point towards the metalicity of NGC2287 being close to the 
solar value, perhaps marginally less. 

The youthful age, the relative proximity and the low line of 
sight reddening of NGC2287 make it a particularly suitable tar- 
get for the study of the IFMR, as was recognised several decades 
ago. Indeed, Romanishin & Angel (1980) undertook a search of the 
cluster using photographic plates and this led to the identification 
of five candidate white dwarf members. Follow-up spectroscopic 
observations obtained with the 3.6m ESO telescope and the Image 
Dissector Scanner confirmed at least three of these objects to be 
white dwarfs, two with Toff^25000K and one with TeS ~13000K 
(Koester & Reimers 1981). It was concluded by these authors that 
the two hotter stars were likely to be cluster members while the 
cooler white dwarf was probably a foreground object. However, as 
the quality of the existing spectral data is rather poor, no detailed 
analyses of the Balmer line profiles in the spectral energy distribu- 
tions of these stars has ever been undertaken. Thus robust estimates 
of the effective temperatures and surface gravities of these objects 
are unavailable to confirm or otherwise these conclusions and to 
allow them to be fully exploited in the context of the IFMR. 



3 OBSERVATIONS AND DATA ANALYSIS 

3.1 FORSI spectroscopy of white dwarf candidate members 
of NGC3532 and NGC2287 

Low resolution, high signal-to-noise optical spectroscopy of the 
nine white dwarf candidate members of the clusters NGC3532 and 
NGC2287 was obtained in service mode with the European South- 
em Observatory (ESO) Very Large Telescope and the Focal Re- 
ducer and low dispersion Spectrograph (FORSI) within the peri- 
ods 2007/04/24-27 and 2007/1 0/06- 11/21. A full description of the 
FORSI instrument may be found on the ESO webpage^ As these 
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targets are comparatively bright, we specified fairly relaxed con- 
straints on the sky conditions and thus the observations were gen- 
erally undertaken in poorer seeing and/or with some cloud present. 
All data were acquired using the 2x2 binning mode of the E2V 
CCD, the 600B+12 grism and a 1.6" slit which gives a nominal res- 
olution of A/AA~500. Flat and arc exposures were obtained within 
a few hours of the acquisition of each of the science frames. 

The CCD data were debiased and flat fielded using the IRAF 
procedure CCDPROC. Cosmic ray hits were removed using the 
routine LACOS SPEC (van Dokkum 2001). Subsequently the spec- 
tra were extracted using the APEXTRACT package and wave- 
length calibrated by comparison with the He+HgCd arc spectra. 
Remaining instrument signature was removed using a spectrum of 
the featureless DC white dwarf WDOOOO+345 obtained with an 
identical set-up during this programme. 

3.2 The model atmosphere calculations 

We have used recent versions of the plane-parallel, hydrostatic, 
non-local thermodynamic equilibrium (non-LTE) atmosphere and 
spectral synthesis codes TLUSTY (v200; Hubeny 1988, Hubeny 
& Lanz 1995) and SYNSPEC (v48; Hubeny, I. and Lanz, T. 2001, 
^http://nova.astro.umd.edu/l to generate a grid of pure-H synthetic 
spectra covering the Toff and surface gravity ranges 14000-35000K 
and log (;=7.25-8.75 respectively. We have employed a model H 
atom incorporating the 8 lowest energy levels and one superlevel 
extending from n=9 to n=80, where the dissolution of the high ly- 
ing levels was treated by means of the occupation probability for- 
malism of Hummer & Mihalas (1988), generalised to the non-LTE 
situation by Hubeny, Hummer & Lanz (1994). The calculations as- 
sumed radiative equilibrium and included the bound-free and free- 
free opacities of the H^ ion and incorporated a full treatment for the 
blanketing effects of HI lines and the Lyman ~a, ~I3 and —7 satel- 
lite opacities as computed by N. Allard (e.g. Allard et al. 2004). 
During the calculation of the model structure the hydrogen line 
broadening was addressed in the following manner: the broaden- 
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ing by heavy perturbers (protons and hydrogen atoms) and elec- 
trons was treated using AUard's data (including the quasi-molecular 
opacity) and an approximate Stark profile (Hubeny, Hummer & 
Lanz 1994) respectively. In the spectral synthesis step detailed pro- 
files for the Balmcr lines were calculated from the Stark broadening 
tables of Lemke (1997). 



Table 2. Coefficients determined for Equation 1, the transformation be- 
tween instrumental magnitudes and Johnson V 



Ko Ki K2 K3 

22.326±0.012 -0.104±0.009 0.213±0.014 -0.138±0.011 



3.3 Determination of the effective temperatures and surface 
gravities 

As is our previous work (e.g. Dobbie et al. 2006a), comparison 
between the models and the data is undertaken using the spectral 
fitting program XSPEC (Shafer et al. 1991). In the present analysis 
all lines from H-/3 to H-8 are included in the fitting process. XSPEC 
works by folding a model through the instrument response before 
comparing the result to the data by means of a x^— statistic. The 
best fit model representation of the data is found by incrementing 
free grid parameters in small steps, linearly interpolating between 
points in the grid, until the value of x'^ is minimised. Errors in the 
TeffS and log g s are calculated by stepping the parameter in ques- 
tion away from its optimum value and redetermining minimum 
until the difference between this and the true minimimi corre- 
sponds to la for a given number of free model parameters (e.g. 
Lampton et al. 1976). The results of our fitting procedure are given 
in Table 1 and shown overplotted on the data in Figure 1. It should 
be noted that the parameter errors quoted here are formal Icr fit er- 
rors and undoubtedly underestimate the true uncertainties. In our 
subsequent analysis we assume more realistic levels of uncertainty 
of 2.3% and 0.07dex in effective temperature and surface gravity 
respectively (e.g. Napiwotzki, Green & Saffer 1999). 
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Figure 2. The derived distance modulii of the 9 white dwarfs included in 
this study. The distance modulus of NGC3532 as estimated by Meynet, 
Mermilliod & Maeder (1993; m-M=8.35), Robichon et al. (1999; m- 
M=8.04) and Kharchenko et al. (2005; m-M=8.61) and the distance modu- 
lus of NGC2287 as estimated by Harris (1993; m-M=9.23), Meynet, Mer- 
milliod & Maeder (1993; m-M=9.15), Kharchenko et al. (2005; m-M=9.30) 
and Sharma et al. (2006; m-M=9.26) are overplotted. The white dwarfs 
which are non-members are clearly distinguished. 



3.4 V band CCD photometry of the white dwarfs 

V band CCD imaging of the nine white dwarfs candidate mem- 
bers of NGC2287 and NGC3532 was obtained on the nights of 
2008/03/06 and 2008/03/07 with the Australia National Univer- 
sity's 40" telescope and the Wide Field Imager (WFI) located at 
Siding Spring Observatory. Conditions on both nights were com- 
paratively good with photometric skies and seeing as measured 
from the images of ~1.5". The WFI consists of a mosaic of 8 
MIT Lincoln Labs 4096 x 2048 pixel CCDs which covers an area 
of 52' x52' in each pointing. Throughout our run, however, CCD6 
was non-functional. All data were reduced using the Cambridge 
Astronomical Survey Unit CCD reduction toolkit (Irwin & Lewis 
2001) to follow standard procedures, namely subtraction of the 
bias, flat-fielding, astrometric calibration and stacking. Aperture 
photometry was performed on the stacked images using a circu- 
lar window with a diameter of 1.5 x the full width half maximum 
of the mean point spread function. The Landolt fields SA98 and 
SA104 (Landolt 1992) were observed a number of times during the 
latter night so that instrumental magnitudes could be transformed 
onto the standard V Johnson system, 

mv = -2.5log{ADU/tea,p) + Kq + KiX + K^^B - V) (1) 

+KzX{B - V) 

where ADU is a measure of the total counts from the source, tcxp 
the exposure time and X the airmass. The coefficients and their 
respective errors were determined to have the values shown in Ta- 
ble 2. The B — V colour of each white dwarfs was estimated from 



the known effective temperature and surface gravity using the syn- 
thetic photometry of Bergeron et al. (1995) as updated by Holberg 
& Bergeron (2006). Our estimates of the V magnitudes of the nine 
white dwarfs are listed in the final column of Table 1. 



4 DISCUSSION 

4.1 Membership status of the nine white dwarfs 

We have used the estimates of the effective temperatures and sur- 
face gravities of the white dwarfs shown in Table 1 and the model 
grids of Bergeron et al. (1995), as revised by Holberg & Bergeron 
(2006), to derive absolute magnitudes (My; see Table 1). Subse- 
quently, we have determined the distance modulii of the nine white 
dwarfs, neglecting extinction which is believed to be low Av ^0. 12 
along the lines of sight to NGC3532 and NGC2287. These are plot- 
ted (solid bars) along with a number of distance estimates available 
in the literature for each of these clusters (dash-dotted lines) in Fig- 
ure 2. It is clear from Figure 2 that NGC3532-6 and NGC3532-8 lie 
beyond NGC3532 and thus are most probably field objects. More- 
over, -3 appears to lie to the foreground of NGC2287 as concluded 
by Koester & Reimers (1981) and is also probably a field white 
dwarf. The remaining six objects have distance modulii which ar- 
gue strongly that they are members of NGC3532 or NGC2287. 
These stars are suitable for placing constraints on the form of the 
IFMR. 
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Figure 3. The locations of the white dwarf members of NGC3532 and NGC2287 in initial mass-final mass space. Data points from a number of other 
populations with metalicities close to the solar value are also shown (see text for details). The theoretical IFMR of Marigo et al. (2007; dot-dashed heavy line), 
the semi-empirical IFMR of Weidemann (2000; heavy dotted line) and the initial mass - core mass at first thermal pulse relation from Karakas, Lattanzio & 
Pols (2002; medium sohd line) are overplotted. The peak in the field white dwarf mass distribution (ztlcr) is represented by the band of grey shading. 



4.2 The IFMR 

4.2.1 Cluster parameters, initial masses and excluded objects 

The masses and cooling times of the six cluster white dwarfs have 
been determined using modem evolutionary tracks supplied by the 
Montreal group (e.g. Fontaine, Brassard & Bergeron 2001). So that 
this work is generally consistent with other recent studies in this 
area (e.g. Dobbie et al. 2006a, Williams & Bolte 2007) we have 
adopted the calculations which include a mixed CO core and thick 
H surface layer structure. The masses and cooling times shown 
in Table 1 have been derived using cubic splines to interpolate 
between points in this grid. We note that these mass determina- 
tions are relatively insensitive to our choice of core composition 
and if instead we had adopted thin H-layer models these estimates 
would be systematically lower by only O.O2M0, which is within 
our present level of precision. 

The lifetime of the progenitor star of each of these six objects 
has been derived by subtracting the estimated cooling time, shown 
in Table 1, from the adopted cluster age, where, for the reasons 
outlined in Section 2, it has been assumed that r=243±40Myrs for 
NGC2287 and r=300±25Myrs for NGC3532. In this calculation, 
we have taken the errors in the cooling times to have magnitudes 
as shown within the brackets of the relevant column of Table 1. 
These are based on more likely levels of uncertainty in our effec- 
tive temperature and surface gravity determinations of 2.3% and 
0.07dex respectively (e.g. Napiwotzki, Green & Saffer 1999). Sub- 



sequently, we have used cubic splines to interpolate between the 
lifetimes calculated for stars of solar compositior0 by Girardi et 
al. (2000) and have constrained the masses of these six progeni- 
tors to the values shown in the final column of Table 1 . We note 
in this context and in the framework of main sequence turn-off 
based cluster age estimates that current eclipsing binary data are 
consistent with aov=0.2 (ie. a moderate level of) convective core 
over-shooting across the broad mass range M~2-3OM0 (Claret 

2007) . However, we caution that the errors we quote here in pro- 
genitor mass are merely approximations since they have not been 
determined through a detailed statistical analysis (e.g. Salaris et al. 

2008) . Nevertheless, their magnitudes should provide a guide to the 
impact of the uncertainties in both the white dwarf parameters and 
the cluster ages which are main sources of error on the final and 
initial masses respectively (Salaris et al. 2008). 

The locations of the six white dwarf members of NGC3532 
and NGC2287 in initial mass-final mass space are shown plotted 
in Figure 3, with data from the extensively studied Sirius system 
(e.g. Barstow et al. 2005) and a number of other clusters which 
have metalicities that are found to be reasonably close to the so- 
lar value (within 30-40%). While Sirius B is not associated with a 
particular star cluster, we include it here since there are relatively 
few objects in the Minit^5M0 regime and the uncertainties on in- 



^ Based on the Anders & Grevasse (1989) definition 
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dividual points here are particularly large. To determine the initial 
and final masses of the white dwarf members of these additional 
populations we have used the same model grids and methodology 
as applied to NGC3532 and NGC2287. For the Pleiades we have 
assumed an age of r=125±25Myrs (e.g. Ferrario et al. 2005) and 
have utilised the white dwarf parameters determined by Dobbie et 
al. (2006a,2006b). In the cases of NGC6819 (Kalirai et al. 2008), 
NGC7789 (Kalirai et al. 2008), the Hyades (Claver et al. 2001), 
Praesepe (Claver et al. 2001, Casewell et al. 2009), NGC6633 
(Williams & Bolte 2007), NGC1039 (Rubin et al. 2008) and Sir- 
ius (Liebert et al. 2005), we have adopted for the cluster age and 
the effective temperatures and surface gravities of the white dwarf 
members, the values listed in the relevant referenced work. 

Although Kalirai et al. (2005) adopted an age of r=650Myrs 
and a substantially subsolar composition for NGC2099, the re- 
sults from two recent spectroscopic studies suggest that the 
cluster has near solar metalicity. Marshall et al. (2005) mea- 
sure [Fe/H]=-l-0.05±0.05 from moderate resolution spectroscopy 
of eight giant members while Hartman et al. (2008) determine 
[M/H]=-l-0.02±0.04 from high resolution spectroscopy of candi- 
date members with reff>4500K. Moreover, during the last fif- 
teen years the bulk of age estimates for NGC2099 obtained us- 
ing theoretical isochrones generated from solar metalicity stellar 
models which include moderate levels of convective core over- 
shooting, have found values within the range r=450-550Myrs e.g. 
r=450Myrs, Mermilliod et al. (1996), r=520Myrs, Kalirai et al. 
(2001), r=450Myrs, Kiss et al. (2001) and r=550Myrs, Hartman 
et al. (2008). Here we have assumed the cluster has solar metal- 
icity and have adopted the mean of the above age determinations, 
r~490±70Myr, where the error bound has been tuned to envelope 
the bulk of these estimates. 

Both Ferrario et al. (2005) and Dobbie et al. (2006a) adopted 
r=158Myrs for NGC2516, a key cluster for constraining the form 
of the top end of the IFMR. This age was drawn from the work of 
Sung & Bessell (2002) and is marginally larger than that adopted 
by Koester & Reimers (1996), r=140Myrs which is from the work 
of Meynet, Mermilliod & Maeder (1993). Kharchenko et al. (2005) 
recently derived r=120Myrs using 2MASS PSC data, but this is 
based on only three cluster stars. Since the work of Ferrario et al. 
(2005) and Dobbie et al. (2006a), a new detailed photometric study 
of the cluster (Lyra et al. 2006), which used the isochrones of Gi- 
rardi et al. (2000), has concluded that r=140Myrs. Looking at all 
these estimates collectively, we conclude that the age of NGC2516 
most probably lies within the range 145±30Myrs. 

We have excluded a number of white dwarf candidate mem- 
bers of these clusters from our subsequent analysis for the follow- 
ing reasons: WD0837-I-218 is more likely to be a field star than 
a member of Praesepe (Casewell et al. 2009), WD0836-I-201 is 
strongly magnetic and may have a substantially different evolu- 
tionary history to that of a typical non-magnetic star (e.g. Wick- 
ramasinghe & Ferrario 2000, Tout et al. 2008), WD0836-I-185, 
NGC6633 LAWDS 4 and 7 may be double-degenerate systems 
as suggested by photometric or radial velocity data, in which case 
close binary interaction could have significantly impacted their evo- 
lution, NGC6633 LAWDS 16 is a DB white dwarf for which de- 
terminations of effective temperature and surface gravity are con- 
siderably less certain (e.g. Kepler et al. 2007), NGC2099 WD 6 
and 21 do not have spectroscopic surface gravity determinations, 
NGC2099 WD 15,16 and 17 are found to be too old at the re- 
vised age of their putative parent cluster, NGC1039 LAWDS 20, 
SI and 9 have proper motions, as listed in the SuperCOSMOS 
Sky Survey database (Hambly et al. 2001), which are ~3.3(7, 
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Figure 4. SuperCOSMOS proper motion measurements for white dwarf 
candidate members of NGC1039 recently identified by Rubin et al. (2008). 
The motions of the three massive white dwarfs (filled triangles) are consis- 
tent with those of the cluster members colour selected by Irwin et al. (2006; 
small filled circles). The white dwaifs with masses close to O.6M0 are also 
shown (filled squai'es). The proper motions of LAWDS 20 and LAWDS SI 
are inconsistent with membership of NGC1039, while the proper motion of 
LAWDS 9 indicates it too is more likely to be a field star. 

~4.1(T ~2.3cr from the mean of their putative parent cluster and 
are thus likely to be field stars (see Figure 4) and NGC7063 
LAWDS 10, with M^0.4Mq (Williams & Bolte 2007), has a 
proper motion as listed in the USNO-B 1.0 catalogue of jia cos 5=- 
2±4 mas yr~^ and /ii=-22±8 mas yr~^ and as measured by 
us from blue Palomar Sky Survey plates (0269, 1951/07/05 and 
SJ04686, 1992/07/26) of cos5=2.1±5.9 mas yr"^ and ^if,=- 
19.0±5.6 mas yr^^, which argues it is more likely to be a field star 
than a cluster member (^a cos 5=1.24±0.41 mas yr^^ and /i5=- 
2.83mas yr"^ Kharchenko et al. 2005). 

4.2.2 A new look at the form of the IFMR 

What is immediately apparent from Figure 3 is that the white 
dwarfs of NGC3532 cannot now substantiate any claim that there 
is significant scatter in the IFMR as could their location in the plots 
of some other recent studies of the IFMR (e.g. Ferrario et al. 2005, 
Catalan et al. 2008). The NGC3532 data points in these investiga- 
tions are based on the older spectroscopy of Reimers & Koester 
(1989) and Koester & Reimers (1993) which is of lower quality 
than the data presented here. At the revised cluster age and met- 
alicity, the locations of the majority of the NGC2099 white dwarfs 
now appear to be entirely consistent with those of the Hyades, Prae- 
sepe, NGC3532 and NGC6633 clusters. There are only six sub- 
stantially deviant data points amongst the forty-nine shown in Fig- 
ure 3. It is noteworthy that it is the white dwarf population of the 

^ The proper motion is listed in the SuperCOSMOS Sky Survey database 
as fia cos (5fa-80±307 mas yr~ ^ and f»-60±302mas yr^ ^ . These huge 
uncertainties may indicate that there was a problem with the measurement. 
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most distant cluster in this sample, NGC2099, where proper mo- 
tions are lacking and distance determinations are most unreliable, 
which presents the highest proportion of outliers, 36% or five out 
of these six. Kalirai et al. (2005) estimate a contamination level of 
~25%, equating to ^^^6 interlopers in their total spectroscopically 
observed sample of 24 objects, so it seems feasible that at least 
some proportion of these deviant points are simply field stars. In- 
deed, the masses of four of these objects reside within the peak in 
the field white dwarf mass distribution, (M=O.565M0, <j=Q.Q%Mq, 
Liebert, Bergeron & Holberg 2005), shown by the grey shaded re- 
gion. Proper motion measurements for white dwarf candidate mem- 
bers of NGC1039 (Rubin et al. 2008) which appeared to occupy this 
part of initial mass-final mass space (Rubin et al. 2008) confirm that 
two stars (NGC1039 LAWDS 20 and SI) and argue strongly that 
a third (NGC1039 LAWDS 9), are simply field objects (see Figure 
4). Moreover, with our improved photometry and spectroscopy of 
the NGC3532 white dwarfs we have also demonstrated here that 
the two candidate members which would otherwise reside in this 
region of initial mass-final mass space are field stars. 

Another possibility that cannot be discounted is that some of 
these low-lying stars have formed through a close binary evolution- 
ary channel, most probably where the post-main sequence phases 
have been terminated prematurely by the loss of envelope mass 
through the formation of a common envelope (Willems & Kolb 
2004). The modelUng of Iben & Tutukov (1985) demonstrates that 
a M=4M0 star which experiences Roche Lobe overflow while as- 
cending the Red Giant Branch can result in the production of a 
Mfinai=O.523M0 white dwarf. It would be informative, although 
difficult in practice, to obtain near-IR photometry to search for 
evidence of the presence of cool low mass companions to these 
NGC2099 white dwarfs once they had been shown to have proper 
motion consistent with the cluster. It has also been suggested a 
number of times that binarity could have played a role in the dis- 
tinctive location of LB5893 above the bulk of stars in initial mass- 
final mass space (e.g. Claver et al. 2001, Casewell et al. 2009). Prae- 
sepe is known to harbour a number of blue straggler stars (e.g. 40 
Cancri and Epsilon Cancri; Ahumada & Lapasset 2007). A propor- 
tion of these likely form as a consequence of mass-transfer between 
and the eventual coalescense of the components of primordial bi- 
naries (e.g. Lombardi et al. 2002). In terms of evolution, these stars 
appear to be retarded with respect to their parent population. Thus 
it seems plausible that the progenitor of LB5893 was a blue strag- 
gler and thus the whole evolution of this star has been somewhat 
delayed with respect to the general Praesepe population. NGC2099 
WDll, which also sits above the bulk of stars in initial mass-final 
mass space, might have a similar evolutionary history. 

We find that the current crop of solar metalicity open clus- 
ter white dwarfs appears to offer no decisive evidence that non- 
magnetic stars which have effectively evolved in isolation, experi- 
ence strong differential mass loss, at least within the initial mass 
range explored here. A number of recent studies have shown that a 
simple linear function, Mfinai=Minit m + c, is a reasonable approx- 
imation to the semi-empirical IFMR. For example, based on data 
from the Sirius binary system and seven open clusters, including 
the metal poor NGC2168, Ferrario et al. (2005), determined best 
fit parameters of m=0.10038±0.00518 and c=0.43443±0.01467 
over the range 2.5M0^Minit;$6.5M0. With the addition of data 
from three more open clusters, Kalirai et al. (2008) obtained pa- 
rameters of m=0.109±0.007 and c=0.394±0.025 over the range 
1.16M0;$Minit^6.5M0. However, there is some evidence within 
the data shown in Figure 3 that the IFMR is somewhat steeper 
in the range 3M0;$Minit;$4M0 than for initial masses immedi- 



ately lower and higher than this. For example, taking the four 
lowest mass Hyades stars and the white dwarfs of NGC7789 
and NGC6819, spaiming the range 1.6M0;$Minit^3.OM0, 
we estimate AMfinai/AMinit=0.0981±0.0301. For the white 
dwarfs lying within the range 3M0S3Minit^4M0, exclud- 
ing the outUers we determine AMfinai/AMinit=0.2279to;j707. 
Over the initial mass regime 3.8M0;$Minit;^5M0 we estimate 
AMfinai/AMinit=0.123llu;o62!, where for the NGC1039 white 
dwarfs alone we find AMHnai/AMinit<0.1269. We note that Fer- 
rario et al. (2005) have demonstrated that population synthesis 
models which adopt an IFMR which has a feature of this nature 
can reproduce the relatively sharp drop observed at the high mass 
end of the main peak in the mass distribution of white dwarfs (e.g. 
Kepler et al. 2007, Liebert, Bergeron & Holberg 2005, Marsh et al. 
1997). 



The trends outlined by the bulk of the data in Figure 3 
bear some resemblance to those of the initial mass-core mass 
at first thermal pulse relation (thin solid line). This theoretical 
track is likewise relatively flatter at Minit^SM© since stars with 
Minit^2.3M0 develop comparable degenerate He cores after core- 
H exhaustion (e.g. Wagenhuber & Groenewegen 1998). It also be- 
comes notably steeper for 3M0^Minit^4M0, as a result of the 
sensitivity of the mass of the H-exhausted core to the initial mass 
(e.g. Becker & Iben 1979). The slope of this relation then decreases 
at Minit^4M0, since here the He-burning shell in the early- AGB 
phase is predicted to be sufficiently potent to power a (second) 
dredge-up event which reduces the mass of the H-depleted core 
(Becker & Iben 1979). The reasons that the IFMR should follow 
rather closely the initial mass-core mass at first thermal pulse rela- 
tion, more especially for Minit^3M0 , have been discussed in depth 
by Weidemann (2000). His estimate of the IFMR is overplotted on 
the data in Figure 3 (heavy dotted line). In brief, during the inter- 
pulse period when the He burning shell is inactive, the H-exhausted 
core increases in mass. However, the re-ignition of the He-shell in 
a flash, drives intershell convection which for a short period ex- 
tends into the CO core, mixing C-rich material up into this zone. 
Subsequently, the H-rich convective envelope extends down into 
the intershell region, dredging processed elements up to the stellar 
surface and effectively reducing the mass of the H-deficient core. 
If dredge-up is particularly efficient, that is, the mass of material 
mixed into the predominantly H-envelope is comparable to the in- 
crease in the mass of the core, A~l, then the core does not grow 
appreciably while the star evolves on the thermally pulsing AGB. 
Detailed modelling of AGB evolution indicates that the maximum 
A value attained during the TP- AGB is a strong function of initial 
mass, with Minit^3M0 reaching large A after only a few thermal 
pulses (e.g. Herwig 2000, Karakas, Lattanzio & Pols 2002). In- 
deed, a recent theoretical IFMR for solar metalicity (Marigo & Gi- 
rardi 2007; dot-dashed heavy line) reproduces the steepening in the 
range 3^Minit;$4M0 and the decrease of the slope at Minit<^4M0 . 
However, the data in Minit^3.3M0 regime appear to sit slightly 
above (a few hundreths of a solar mass) both this and the initial 
mass-core mass at first thermal pulse relation which may indicate 
that third dredge-up may not be quite as efficient here as assumed 
in the Padova models. Nevertheless, the similarities between the 
forms of the theoretical relations and the trends delineated by the 
bulk of white dwarfs from solar metalicity open clusters lend some 
assurance to the results of modem stellar evolutionary calculations. 
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5 SUMMARY 

We have obtained high signal-to-noise low resolution optical spec- 
troscopy of the nine candidate white dwarfs members of NGC3532 
and NGC2287 with FORS 1 and the VLT. The analysis of these data 
and of new V band photometry indicate that only six of these ob- 
jects are probably members of the clusters. These six objects, in 
particular the four members of NGC3532, do not substantiate any 
claim that there is substantial scatter in the IFMR. While a simple 
linear fit to these data could still be deemed acceptable, there are 
now clear hints that the IFMR is steeper in the initial mass range 
3M0^Minit^4M0 than at progenitor masses immediately lower 
and higher than this. This result is consistent with the predictions 
of stellar evolutionary models. Moreover, it can help explain the 
relatively sharp drop in the number density of white dwarfs on the 
high mass side of the main peak in the white dwarf mass distri- 
bution. Unfortunately, the IFMR remains rather poorly constrained 
at Minit^5MQ, where there is particular interest in its form. Ad- 
ditional white dwarfs and improved spectroscopy on exisiting data 
points are urgently required in this progenitor mass range. 
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